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Abstract We present a model to estimate the synchrotron 
radio emission generated in microquasar (MQ) jets due to 
secondary pairs created via decay of charged pions produced 
in proton-proton collisions between stellar wind ions and jet 
relativistic protons. Signatures of electrons/positrons are ob- 
tained from consistent particle energy distributions that take 
into account energy losses due to synchrotron and inverse 
Compton (IC) processes, as well as adiabatic expansion. The 
space parameter for the model is explored and the corre- 
sponding spectral energy distributions (SEDs) are presented. 
We conclude that secondary leptonic emission represents a 
significant though hardly dominant contribution to the to- 
tal radio emission in MQs, with observational consequences 
that can be used to test some still unknown processes oc- 
curring in these objects as well as the nature of the matter 
outflowing in their jets. 
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1 Introduction 



X-ray binary systems (XRBs) are composed by either a stel- 
lar mass black hole or a neutron star, and a normal (non 
degenerated) star which supplies matter to the compact ob- 
ject through the formation of an accretion disk. Some 260 
XRBs are known up to now [1] probably corresponding to 
an underlying population of some tens of thousands of com- 
pact objects in our Galaxy. A few of these sources present 
also non-thermal radio emission, hence evidencing the ex- 
istence of mechanism(s) capable of injecting and/or accel- 
erating large numbers of relativistic particles. Some radio 
emitting X-ray binary systems (REXBs) have been observed 
showing ejection of material at relativistic velocities and to 
display jets like those seen in quasars and active galactic nu- 
clei but at ~ 10~ 6 times shorter scales. This analogy is the 
reason for calling them microquasars (MQs) [2] and making 
them some of the most interesting objects for astrophysics. 
Furthermore, attention on these objects has grown since the 
proposal of Paredes et al. (2000) [3| of MQs as counter- 
parts of some of the unidentified gamma-ray sources of the 
EGRET catalog [4] and hence pointing them as plausible 
high energy emitters. A strong confirmation of this associa- 
tion has come from the detections of the MQs LS 5039 and 
LS I +61 303 at Tev energies using respectively the ground- 
based Cherenkov telescopes HESS [5 ] and MAGIC |6), giv- 
ing support and empowering at the same time a number of 
previous detailed studies centered on the mechanisms op- 
erating in these sources in order to explain the gamma ray 
domain (see, e.g., [7 | and [8 |). Moreover, a jet origin of the 
emission from MQs has been suggested from the observa- 
tion of syncrothron emission of relativistic electrons/positrons 
extending from the radio all the way into the X-ray regime. 
In this sense jet-like models have focused on different ap- 
proaches regarding the particle origin that could generate 
the required emission properties in a consistent way. Some 
of them consider leptons directly injected at the base of the 
jet and, in extending outwards, Compton-interaction with 
external/self-created photon fields produces high energy ra- 
diation. Other models deal with an hadronic origin of the 
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high energy emission, through proton-proton interactions and 
pion decay producing gamma rays and leaving the result- 
ing co-generated leptons as low energy emitters. The present 
work refers to the later procedure, focusing on the modelisa- 
tion of the secondary leptonic synchrotron emission in order 
to constrain the characterization of MQ jets. An outline of 
the model is given in the next, followed by the results show- 
ing the SEDs and lightcurves under different parameter as- 
sumptions and the conclusions that can be extracted from 
them. 



2 Model description 

We have taken a binary system formed by a black hole and a 
high-mass early-type star which feeds the accretion trans- 
fer of mass onto the compact object while developing an 
accretion disk. Part of the accretion power is converted to 
kinetic and magnetic energy of the accretion flow under the 
effects of the compact object potential well and a relativistic 
e ± -p plasma is ejected in a direction taken to be perpendic- 
ular to the plane defined by the accretion disk. This plane 
is the same as the orbital one, even if this condition can be 
relaxed to allow a more general situation but for simplic- 
ity we assume coplanarity in our model. Jet energetics is 
assumed to be dominated by accretion, and further energy 
sources like compact object rotation have been neglected at 
this stage. The jet will also contain a magnetic field Bj et as- 
sociated with the plasma. We assume that the matter kinetic 
luminosity is higher than the magnetic luminosity (or total 
magnetic energy crossing the jet section per time unit) in the 
jet regions we are concerned with, although the magnetic 
field can be still significant once the jet is formed, since the 
ejection mechanism have likely a magneto-hydrodynamical 
origin. 

In this scenario we deal with radiative processes that take 
place in the jet and can produce significant emission in the 
radio spectrum. Although the contribution of protons could 
still be significant from the radiative point of view, we will 
focus on the leptonic component only. The reader is related 
to other works (see, e.g., |9) and flU) for treatments on pri- 
mary hadrons. 



Table 1 Parameter values used throughout the model. 
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Black hole mass 
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Injection point 
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5QR g 


Initial radius 


Ro 




Radius of the companion star 


K* 


irn 

15K e 


Orbital radius 


a 




Luninosity companion star 
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1.6 x 10 39 erg/s 


Mass loss rate 


M t 


S-ltr^Moyr- 1 


let's Lorentz factor 


r 


1.02 


let kinetic luminosity 


Qj 


10 36 erg/s 


Proton kinetic Luminosity 


Qp 


10 35 erg/s 


Minimum proton energy 


rpinin 


2 GeV 


Maximum proton energy 


T?max 


100 TeV 


Leptonic spectral index 


p 


2.2 


Wind-jet penetration factor 


f P 


0.1 


Magnetic field at zo 


Bo 


10 2 ,10 3 , 10 4 



jet region penetrates into it. When charged and neutral pions 
are created, the first will decay to muons and subsequently 
to electrons and positrons, while the second decays to high 
energy photons through 
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For an injection proton spectrum as the one given above, the 
pion spectrum (in the jet reference frame) will be a power 
law and the electron/positron distribution will also follow a 
power law [10] with a differential pair injection rate given 
by J e {E e ) = K e E~ p where we use a value P = 2.2. 



2.2 Secondary evolution 

This injected electron-positron population suffers subsequent 
radiative cooling due to synchrotron and inverse Compton 
losses as well as adiabatic expansion given by 
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2.1 Secondary generation 

Secondary leptons are generated in p-p interactions of hot 
protons in the jet with the companion star wind that extents 
out isotropically producing charged and neutral pions (7T ± , 
7T )through the reaction channel p +p — > p + p + ^ n o7t° + given by U star = 



dE e ± ' 




dt 






sync 


dE e± ~ 




dt 


IC 


dE ± " 




dt 




exp 



% Jl ±(n + +1t ), where t, K ± ~ 2(E p jGeV) x l A is the charged 



pion multiplicity. The relativistic injection proton spectrum 
is a power law N p (E p ) — K p E~ a where E""" < E 



respectively. Here U s tar is the companion star's photon field 
where d is the distance from the star 
10 39 erg/s, and V exp = dR(z)/dt is 



where the constant K p can be found normalisation to the to- 
tal power that goes to protons Q p . The corresponding proton 
flux will be given by J p (E p ) = (c/4tz)K p {z /z) 2 E p a . We 
assume that a fraction of 1/10 of the matter that crosses the 
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and in our case L+ = 1 .6 ±u ^ mm , exp 
the lateral expansion velocity of the jet. Fresh electrons and 
positrons that have been injected at a certain point into the jet 
suffer the different energy losses modifying in this way the 
spectral distribution of particles at each height z. Moreover, 
different evolution stages sum at each height since injection 
of fresh particles occurs all along the jet and therefore a mix 
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Fig. 1 The strength of the magnetic field B is treated as a free pa- 
rameter in our model. Synchrotron loses will be enhanced when taking 
higher values for B but, since injection occurs all along the jet inside the 
bynary system (until a distance roughly given by z — R or b ~ 10 13 cm), 
the enhancement in the final emission will preval over losses. Here we 
show the SEDs for three different magnetic fields (B); values indicated 
are taken at the base of the jet. 



of multiple evolved population of electron/positron distribu- 
tion arises. To compute it in a consistent way, one has to 
assume the continuity equation, 

N e ± (E e ± , z)dE e ± = N Qe ± (E 0e ± , zo)dE ^ e ± 



where Nn 



■ Kq e ±E„ p ± is the initial energy spec- 



trum for the injected particle density. Solving the equations 
for the evolution of the particle energy along the jet axis, 
one can find the spectral distribution at each slice, compute 
de differential luminosity at each height and finally integrate 
to find the total luminosity. 



2.3 Secondary emission 

Secondary electrons and positrons will radiate through syn- 
chrotron process since we assume the presence of a mag- 
netic field wich is tangled to the plasma and oriented ran- 
domly in direction. The resulting emission will be isotropic 
in the jet reference frame; to compute the expected lumi- 
nosity, we have used the formulae presented in 1111 . where 
expressions for the specific emission and absorption coeffi- 
cients are given by 



£ v = C 3 B £N e ±(E e ±)F(x)dE e 



7J> C 3 B Jo ' E e±JE 



F(x)dE e 



with the function F(x) given by F(x) = x K 5 ^(z)dz be- 
ing K 5 p the modified second kind Bessel functions, x = ^- 
where V c = T ^ T ^E 2 ± « 6.27 • 10 18 £ 2 ± Hz is the critical fre- 

1.87 ■ 10 -23 (in CGS units). 
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Fig. 2 Light curves of emission at 1 and 5 GHz bands at growing dis- 
tances along the jet for various magnetic fields B wich values are ref- 
ered to the base of the jet at z = 50R g . Not only the overall emission 
but also the peak of the lightcurves moves to higher distances when in- 
creasing B, as well as the difference between the locations of the peaks 
for the two frequencies. 
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Fig. 3 Variations of the SEDs due to different observing angles 6 be- 
tween the approaching jet axis and the observer line of sight.As far as 
the jet material travels at higher velocities along the jet, the boosting 
will get increasingly greater while observing at smaller angles. The 
SEDs showed here correspond to 6 ranging from to j , where we 



have used a Lorenz factor F ■ 
tivistic jet with vj ~ 0.2c 



1.02 corresponding to a mildly rela- 



Then, the specific differential luminosity at a certain heigh z 
in the jet can be expressed as 

^=2*R( Z )£(l-«-'A) 

where I; is the typical size of the synchrotron emitting plasma 
region. Now, integrating over the entire jet lenght we find the 
spectral energy distribution, 



quency, and the constant C3 



v jz*™ §2dL0l dz 
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where 8 = [r (1 — j8 cos 9 t, s )] 1 is the Doppler boosting fac- 
tor. 



3 Model results and conclusions 

SEDs are obtained for different magnetic field values, elec- 
tron/positron spectral indices and spatially distributed disks. 
We have estimated also the expected emission along the jet 
at 1 and 5 GHz. Leptons are injected in the context of hadronic 
secondaries generation within a detailed model that takes 
into account in a consistent way particle injection mecha- 
nisms and cooling due to radiation processes and adiabatic 
expansion. The luminosities obtained are slightly lower than 
in the models based on primary leptons injection, and must 
be considered complementary to them. However, we note 
that within our model there is no requirements of accelera- 
tion processes along the jet to obtain the final emission re- 
sults. Such acceleration processes are still not well under- 
stood, although. They could come from diffusive shock ac- 
celeration along the jet when fresh ejecta interact with previ- 
ous blobs of plasma already outflowing at lower velocities. 
Other scenarios assume a continuous energy transfer mech- 
anism from the magnetic field to the matter content of the 
jet in such a way that the resulting parsec-scale radio emis- 
sion can be explained. The fact of studying alternative mod- 
els were particles are directly injected until a certain height 
along the jet can constrain the amount of acceleration re- 
quired and contribute to the understanding of the physical 
mechanisms that can lead to such processes. 

Signatures at different distances along the jet and spe- 
cific spectral features detectable for reasonable parameter 
values treated in our numerical simulations have the poten- 
tial to be an important clue for determining the matter con- 
tent of jets. In particluar, highly resolved observations at 1 
and 5 GHz could determine if leptons are present at heights 
1Q12-13 cm at tne ec [g e f the binary system typical region 
where wind matter from the companion is still significant. If 
electrons/positrons still show high energies due to a recent 
injection from hadronic interactions at these parts of the jet, 
it could be a signature of secondary generation without the 
necessity of invoking additional acceleration processes. 
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